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The anharmonic decay rates of atomic vibrations in amorphous silicon (a-Si) and paracrystalline 
silicon (p-Si), containing small crystalline grains embedded in a disordered matrix, are calculated 
using realistic structural models. The models are 1000-atom four-coordinated networks relaxed 
to a local minimum of the Stillinger- Weber interatomic potential. The vibrational decay rates 
are calculated numerically by perturbation theory, taking into account cubic anharmonicity as the 
perturbation. The vibrational lifetimes for a-Si are found to be on picosecond time scales, in 
agreement with the previous perturbative and classical molecular dynamics calculations on a 216- 
atom model. The calculated decay rates for p-Si are similar to those of a-Si. No modes in p-Si reside 
entirely on the crystalline cluster, decoupled from the amorphous matrix. The localized modes with 
the largest (up to 59%) weight on the cluster decay primarily to two diffusons. The numerical results 
are discussed in relation to a recent suggestion by van der Voort et al. [Phys. Rev. B 62, 8072 
(2000)] that long vibrational relaxation inferred experimentally may be due to possible crystalline 
nanostructures in some types of a-Si. 

PACS numbers: 63.50.+x,65.60.+a 



I. INTRODUCTION 

The pioneering experiments by Dijkhuis and 
coworkersi^*^ explored transient dynamics of ex- 
cited vibrational modes in a topologically disordered 
material-hydrogenated amorphous silicon. In these 
experiments nonequilibrium vibrational states were gen- 
erated during relaxation and recombination of optically 
excited electrons, and monitored with a probe laser 
(anti-Stokes Raman spectroscopy) for transient behav- 
ior. The experimental results are surprising: Scholtcn 
et ali£ found that at low temperatures (2 K) and for 
vibrational frequencies greater than 10 meV (maximum 
frequency in a-Si is about 70 meV) vibrations decay 
on time scales of tens of nanoseconds. Furthermore, 
the higher the vibrational frequency, the slower is the 
decay rate. In contrast, phonons in crystalline silicon 
decay on time scales of tens of picoseconds^ and the 
decay rates increase with increasing frequency. The 
results of Scholten et al. were further confirmed by 
van der Voort et al.^ who suggested that the long 
lifetimes are due to the microstructure of amorphous 
silicon. This suggestion was tested by van der Voort 
et al^ by measuring the vibrational decay rates of a 
mixed amorphous-nanocrystalline silicon, which was 
an amorphous hydrogenated silicon with a sizable 
fraction of nanocrystallites (with the diameter of 1 — 5 
nm). Even the mixed sample displayed nanosecond 
vibrational lifetimes, although the lifetimes appeared to 
decrease with increasing frequency. A hypothesis was 
put forward^ that the measured types of amorphous 



silicon contain nanoscale regions with correlated (if 
not ordered) atoms, which, through enhanced size 
quantization and localization of vibrational frequencies, 
inhibit anharmonic decay. 

These experimental results are at odds with the known 
theories of anharmonic vibrational decay in disordered 
materialsAL^ In their so called "fracton" model, Alexan- 
der et ali^ assumed that the majority of vibrational states 
in disordered systems are localized. This seemed to ex- 
plain the above experimental findings since the anhar- 
monic decay could be drastically reduced by the ex- 
tremely small likelihood of the overlap between three lo- 
calized modesii2iii That the small probability of the over- 
lap between three localized modes inhibits vibrational 
decay was disputed by Fabian and Allen^ who put for- 
ward a probabilistic scaling argument that the interac- 
tion between three localized modes would in fact be cru- 
cial for the anharmonic decay and cannot be neglected. 
Fabian later demonstrated^ the scaling argument on a 
one dimensional anharmonic chain with random spring 
constants, and similar conclusions were reached recently 
by Lcitncr in a study of heat flow in a one dimen- 
sional glassi^ and vibrational energy transfer in helices of 
myoglobin^ Thus the fracton model, even if true in its 
premise of localization of the majority of the vibrational 
modes, does not explain the experiment. We note, how- 
ever, that even the premise of the model is questionable, 
as it is in sharp contrast to what is found in finite-size 
realistic models of glasses, which normally exhibit local- 
ization only in a small part of the spectrum. 

Numerical calculations of vibrational decay in glasses 
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have been performed both by evaluating a perturbation 
formula^ and by classical molecular dynamicsi^ Pertur- 
bation theory was applied to the problem of anharmonic 
decay in glasses by Fabian and Allen& who computed 
the decay rates for a 216 atom model of amorphous Si. 
The decay rates were found to be fractions of meV (that 
is, lifetimes are picoseconds), in general increasing with 
increasing frequency. The anharmonic lifetimes of local- 
ized modes were similar to those of the extended modes, 
even in the case of a model alloy Si^Gei-^, where lo- 
calized modes span more than a half of the spectrum 
and the overlap between localized states becomes impor- 
tant£ Bickham and Feldmani reported vibrational decay 
rates for selected modes of 216 and 4096 atom models of 
amorphous Si, using classical molecular dynamics. Their 
results agree with the perturbative calculation, though 
the computed decay rates are somewhat greater due to 
the fact that molecular dynamics takes into account all 
the anharmonic interaction, while the perturbative cal- 
culation in Ref. |6| only cubic anharmonicity. In the cal- 
culation of Bickham and Feldman, a chosen vibrational 
mode was given a greater than average kinetic energy and 
was allowed to equilibrate while keeping the overall tem- 
perature constant. From the decay of the kinetic energy 
in time, the mode decay rate was obtained. While the 
advantage of molecular dynamics over perturbation the- 
ory in calculating vibrational decay rates is that the full 
anharmonic interaction is considered, the disadvantage is 
that the classical dynamics does not capture accurately 
the low temperature decay rates (for example, the rates 
computed by classical molecular dynamics vanish at zero 
temperature^ while in reality they are finite due to quan- 
tum effects^). 

The purpose of this paper is twofold: (i) to extend 
the previous numerical studies of perturbative anhar- 
monic decay in homogeneous amorphous silicon (a-Si) 
to a larger system, and (ii) to calculate vibrational de- 
cay rates for a model of amorphous silicon — paracrys- 
tallinc silicon (p-Si) — that includes nanocrystallites. 
The larger system is a 1000-atom model of a-Si, prepared 
similarly to the previously used 216-atom models The 
calculated decay rates display smaller statistical fluctu- 
ations and agree, on average, with those of the smaller 
model. Also, the larger system has a lower minimum 
frequency, allowing to calculate decay rates at the fre- 
quencies previously unattainable. Studying paracrys- 
tallinc silicon, a material where small crystalline grains 
are embedded in a disordered matrix, allows us to test 
the hypothesis of van der Voort^ regarding the struc- 
tural origin of the anomalous long vibrational lifetimes 
in a mixed amorphous-nanocrystalline Si system. In our 
calculations we have used a 1000-atom (86 out of which 
belong to a single crystalline grain) model created by 
Nakhmanson et. aliU to simulate medium-range order in 
amorphous silicon. We should point out, that, although 
providing a more realistic subject for the verification of 
van der Voort's hypothesis than "regular" models for a- 
Si, this simple model is neither an exact structural match 



to nanocrystalline Si sample of Ref. S (24% crystalline 
fraction and 4.5 nm average grain diameter, versus « 10% 
and 1 nm in the model) nor can it account for various 
other topological defects present in real material. Still, if 
van der Voort's supposition were correct, we would ob- 
serve inhibited decay rates of the modes which would be 
predominantly localized on the crystalline cluster in the 
model. However, we do not find any modes localized ex- 
clusively on the cluster: one of the most cluster-localized 
modes has only 59% weight on the cluster, and is there- 
fore well coupled to the disordered matrix. It is not sur- 
prising that such modes have decay rates similar to other 
localized states. 

We remark that the names amorphous and paracrys- 
tallinc in reference to our models are a matter of termino- 
logical convenience rather than an attempt in classifying 
real materials. We refer to a-Si as describing a homoge- 
neous, continuous random network of silicon atoms, while 
p-Si models are such networks filled with crystalline clus- 
ters. Real materials — which are normally termed amor- 
phous silicon — are likely of the p-Si type, containing 
nanoscale crystallites with a distribution of sizes^ 

In the following we first introduce the structural mod- 
els of a-Si and p-Si and their harmonic vibrational prop- 
erties, then discuss the perturbative calculation of anhar- 
monic decay rates and present the results for the 1000- 
atom models of a-Si and p-Si. Finally, we discuss our 
results with respect to the experiment. 



II. MODELS 

Both models employed in our studies were created 
with similar techniques: the homogeneous model for a- 
Si was constructed using the WWW methodic and the 
paracrystalline model with a variation" of the Barkema- 
Mousseau method^ (For a recent review of modeling 
continuous random networks see Ref. 0.) The for- 
mer model was also studied^ for its harmonic proper- 
ties within the framework of the Stillinger- Weber (SW) 
potential2i prior to the present work. Apart from the 
computational efficiency issues, the major difference be- 
tween the methods of WWW and Barkema-Mousseau is 
the starting configuration used for the model construc- 
tion: crystalline silicon (c-Si) is used in the former ap- 
proach and a random close packed configuration in the 
latter. Since neither of the starting models of this work 
employ the SW potential, both models had to be relaxed 
with respect to SW prior to the decay times calculation. 
This was done in an earlier work for the homogeneous 
model and in the present work for the paracrystalline 
model through the use of a molecular dynamics quench, 
where both temperature and virial pressure are set to 
zero in this stage of model preparation. The densities for 
the two SW relaxed models are thus found to be slightly 
(three to four percent) less than that of the density of 
c-Si, with the density of the paracrystalline model being 
slightly (two percent) higher than that of the homoge- 




FIG. 1: The 1000-atom model of p-Si. Atoms belonging to 
the crystalline grain are shown in grey; those in white form 
the amorphous matrix. 

neous model. Changes in the atomic positions resulting 
from the SW relaxation were found to be quite small 
and the view of the SW relaxed model of p-Si shown in 
Fig. ^ is very close to the unrelaxed structure. In gen- 
eral it is known that the SW potential produces relaxed 
structures that have two to three percent five-fold coordi- 
nated atoms, even if the starting structures were perfectly 
four-fold coordinated (which is the case for the paracrys- 
talline model). At variance with the optical properties 
of the models, this deviation from the perfect four-fold 
coordination does not noticeably alter their vibrational 
properties and will not be discussed here in more detail. 



III. HARMONIC VIBRATIONS 

In the harmonic approximation vibrational eigenfre- 
quencies ui(i) and eigenvectors are computed by diag- 
onalizing the corresponding dynamical matrix (through- 
out the paper symbols j, fc, and I will represent vibra- 
tional modes, while a, 6, and c atoms). The results of 
numerical calculations from various g r0 ups^2i2&2£i2£iSii2£ 
indicate that vibrational eigenstates in glasses belong to 
one of four groupa&S&Si: propagons, resonant modes, 
diffusons, and locons. Propagons are essentially sound 
waves scattered by structural disorder. Resonances 
(sometimes called quasilocalized modes) are modes tem- 
porarily trapped in topological defects (in Si models these 
are the groups of undercoordinated atoms) . Spatially, a 
resonant mode has a large vibrational amplitude at the 
defect, but although being rather weak, the mode does 



not decay exponentially (as a truly localized mode would) 
outside of the defect region. In models of a-Si resonances 
often have frequencies below the lowest frequency (which 
is not zero due to the finite size) of the corresponding 
models of crystalline silicon. The low frequency spec- 
tral region of propagons and resonances is a subject of 
great interest in relation to the so called boson peak^S 
which was also investigated in models of a-Si. 31,32 Diffu- 
sons make up the majority of the spectrum. They have 
frequencies above the Ioffe-Regel limit (in principle, as 
suggested in Ref. 13.1 there can be two such limits, one 
for transverse and one for longitudinal acoustic modes), 
so that they cannot be characterized by a wave vector, 
but nevertheless are extended modes and contribute to 
thermal conduction within Kubo theory as shown by the 
numerical calculations of Refs. l2fll34L Finally, locons are 
localized modes in the sense of strong (Anderson) local- 
ization. In most studied models locons form only a small, 
high-frequency part of the spectrum. In our models of a- 
Si many of the most localized locons reside in regions 
where there are five-fold coordinated atoms. 

Experimentally the character of the atomic vibrations 
in glasses has been studied by inelastic x-ray scattering 
in various glassy systems ^^SiSLSi^Si The general con- 
census seems to be that sound waves cease to propagate 
early in the spectrum (for example, at 9 meV in den- 
sified silica^ and at 5 meV in amorphous selenium 3 ^). 
This confirms the existence of the Ioffe-Regel crossover in 
glasses and the propagon picture. The issue of the nature 
of the vibrational modes above the crossover frequency 
(where the models predict diffusons) has not been re- 
solved experimentally yet. The recent experimental and 
theoretical progress is reviewed in Ref. [2^. 

In Fig. |21 we plot the calculated vibrational density of 
states (VDOS) for the models of a-Si and p-Si. Both 
curves look very similar, which is in agreement with 
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FIG. 3: Inverse participation ratio l/p of the vibrational 
states in the 1000-atom models of a-Si (shaded area) and p- 
Si (line without shading). The modes with the frequencies 
above 72 meV (the vertical line) can be considered localized. 
Quasilocalization occurs at low frequencies (resonant modes) 
and around 30 meV which corresponds to the band edges. 
The inset is a detailed view of the low- frequency region. 
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FIG. 4: Weight of the modes at the crystalline cluster as a 
function of mode frequency in the 1000-atom model of p-Si. 
Plotted is the square of the atom displacement summed over 
the atoms forming the cluster. The horizontal line shows a 
weight of 0.086 (8.6%) indicating an unbiased displacement 
pattern. The histogram is VDOS for the cluster (see text). 



the VDOS calculation of Ref. ^| made with a modified 
version of the SW potential^ The calculated spectrum 
agrees rather well with the experimental one^ except 
that the calculation overestimates the highest frequen- 
cies by about 15%. This is a known artifact of the SW 
potential. 

Localization properties of the modes can be judged 
from the participation ratio p(J), which indicates how 
many atoms "participate" in vibrational eigenmodes j. 
Inverse participation ratio 1 jp for a-Si and p-Si, as a func- 
tion of mode frequency is shown in Fig. [3] The majority 
of vibrations in both models is delocalized, with the lo- 
calization transition taking place at around 72 meV (the 
mobility edge). The modes around 30 meV and some 
modes below 10 meV appear to be localized too. The 
latter are resonant modes. The extended modes below 
about 15 meV are propagons, while all the rest are dif- 
fusons (with possibly some longitudinal propagons left 33 
at small frequencies). The localization character in both 
a-Si and p-Si models is similar. The presence of the crys- 
talline cluster does not lead to additional localized modes 
elsewhere in the spectrum. 

In order to understand what fraction of each mode re- 
sides on the crystalline cluster in the p-Si model, we com- 
pute the weight each mode has on the cluster (that is, we 
sum \e J a \ 2 for each j over all atoms a from the cluster). 
The result is shown in Fig. 01 together with a histogram 
of VDOS of the cluster calculated by solving the dynam- 
ical equations for the cluster atoms with the surround- 
ing atoms held fixed. An unbiased mode has a weight 
of 0.086 (8.6%), corresponding to the percentage of the 
atoms making up the cluster. For all the modes below the 
mobility edge the weight fluctuates around 0.086, show- 
ing no special affinity for the cluster. Localized modes, 



as would be expected from their idiosyncratic character, 
can be localized (fully or partially) on, through, or off the 
cluster. None of the modes is localized fully on the clus- 
ter. There are 4 locons with the weight on the cluster of 
30% or greater, the maximum weight being that of 59% 
for a mode with frequency u> — 73.05 meV and partici- 
pation ratio p — 13. The second most localized mode on 
the cluster has the frequency of 72.69 meV, the weight 
of 55% and p = 12. The third and fourth modes are 
more delocalized, having frequency (weighty) lj = 70.67 
(31%, 160) and w = 71.12 (30%, 117), respectively. All 
four modes lie in the mobility edge region. In addition to 
these, there are modes with frequencies around 30 meV 
which have enhanced affinity for the cluster (see Fig. 0J. 
The weight of these modes at the crystalline cluster does 
not exceed 30%, but six modes have the weight between 
20 and 30%. 

Harmonic vibrations in a-Si explain well many ob- 
served thermodynamic properties^ of the material, as 
well as kinetics such as heat flowi22iM Anharmonicity 
does not directly affect heat flow in dielectric glasses, but 
is very important in relaxing the perturbed vibrations to 
maintain local equilibrium (temperature gradient, to be 
specific). More directly, anharmonicity affects thermal 
expansion and sound attenuation. The 1000-atom model 
of a-Si was employed to demonstrate the importance of 
thermal vibrations in both of these phenomena^S^ 3 . It 
was found that anharmonicity is rather weak in a-Si, al- 
though somewhat stronger than in c-Si, primarily due to 
strong anharmonicity of resonant modes. Indeed, res- 
onant modes show giant Gruneisen parameters in the 
model, strongly enhancing the effects of anharmonicity, 
although still within the limits of perturbation theory 
based on cubic anharmonicity. 
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FIG. 5: Calculated decay rate of the mode (a propagon) with 
frequency ui = 12.49 meV in a-Si as a function of temper- 
ature for different widths w of the rectangle function 6 w {oj) 
representing the delta function in Eq. Q The curves are for w 
equal 0.12 (dashed line) and 0.2, 0.4, 2.4 meV (solid lines), 
the order being not mirrored in the magnitude of the curves. 
The greatest decay rate is for to = 0.2 meV, while the lowest 
for w = 1.4 meV. The curve representing w = 1 meV chosen 
in the calculation is in the middle of the bunch. 



IV. VIBRATIONAL LIFETIMES 

Using cubic anharmonicity as the small perturbation to 
the harmonic Hamiltonian, anharmonic decay rate 2T(j) 
of mode j can be obtained from the formula^^ 



2r(j) - 
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n(k) + n(l)} 5 [w(i) - w{k) - oj(1)} 
-n(l)]6[u(j)+u(k)-u>(l)]). (1) 



Here u>(j) is the frequency of mode j, n(j) is the mode 
occupation number given by n(j) = {exp[huj(j) / ksT] — 
with T denoting temperature, and V(j, k, I) is the 
matrix element of the cubic anharmonicity of the inter- 
atomic potential V in the harmonic representation: 
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Greek symbols a, f3, and 7 stand for the cartesian co- 
ordinates of both the atomic displacements u from the 
equilibrium positions, and of the normalized vibrational 
eigenvectors e. The atomic masses are denoted as to. 
Anharmonic vibrational lifetimes are the inverse of the 
rates: 



r(j) = 1/2TU). 



(3) 



In this paper we present decay rates in the units of meV. 
For conversion into lifetimes, a decay rate of 1 meV is 
equivalent to a lifetime of about 0.7 ps. 
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FIG. 6: Calculated decay rate of the mode with frequency 
uj — 12.49 meV in a-Si as a function of w at 10, 100, and 300 
K. 



In Eq. ^ the term with the temperature factor 1 + 
n(k) + n(l) corresponds to the decay j — > k + I, while 
the term with n(k) — nil) represents the "difference" de- 
cay j + k — * /. Energy conservation is ensured by the 
delta functions. At low (down to zero) temperatures the 
first term in Eq. ^ dominates, giving rise to a constant 
2r, while both terms are generally equally important at 
large temperatures, where T ~ T, In crystals V(j, k, I) 
vanishes unless the modes' momentum is conserved in 
the decay process. In glasses, where lattice momentum 
itself is not a valid concept (except for propagons and 
resonances), all the modes k and 1 from the spectrum 
contribute to V(J, k, I) for a given j. 

Crucial for determining decay rates perturbatively 
from a finite-size model is the delta-function regular- 
ization. We approximate 8(u>) ~ 9 w (u>), where 9 w (lo) 
is a rectangle of width w and height 1/w centered at 
uj = 0. In our calculations with 1000 atoms we choose 
w = 1 meV, which fits about 40 modes in the rectangle. 
The choice of w needs to be a compromise between good 
statistics and computer power. The statistics is deter- 
mined by both the number and "similarity" of the modes 
in a rectangle. If w is too large, the rectangle function 
will sample modes with distinct characteristics, not rep- 
resenting faithfully the modes of the chosen frequency. 
This problem is likely to be absent for diffusons, which 
do not differ much on small s pect ral scales due to the 
absence of degeneracy (cf. Ref. l42j) . but may be relevant 
for propagons (which are mixed with resonant modes) 
and locons (which are idiosyncratic). Fortunately, the 
averaging, first within the rectangle and second, over the 
whole spectrum (see Eq. makes the decay rates quite 
insensitive to the choice of w, for a reasonable interval 
of values. In the earlier calculation^ w was chosen to be 
0.4 meV for a 216-atom model, fitting about 4 modes in 
the rectangle. As we will see from the comparison of the 
two calculations in the following section, this choice was 
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FIG. 7: Calculated decay rates versus frequency for the 1000- 
atom model at 10 K (thick line). For comparison the rates 
of the 216-atom model from Ref. |g (not smoothed) are also 
shown (shaded area). 
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FIG. 8: Calculated decay rates versus frequency for the 1000- 
atom model at 300 K. For comparison the rates of the 216- 
atom model from Ref. (not smoothed) are also shown 
(shaded area). 



already good enough, although it may have contributed 
somewhat to the statistical noise, especially at low tem- 
peratures and small frequencies. To illustrate the effect w 
has on the decay rates, we show in Fig. [5] the calculated 
rates of the mode with oj — 12.49 meV in a-Si, as a func- 
tion of temperature, for selected w, ranging from 0.12 
meV (corresponding to about 4 modes per rectangle) to 
2.4 meV (80 modes/rectangle). Except for the smallest 
w, the results are grouped together with the dispersion 
of less than 10% above 100 K. The greatest dispersion is 
at the lowest temperatures, where it reaches 25%. (The 
low temperature properties of the model do not describe 
well the real a-Si structure, because of the existence of 
the minimum frequency of 4 meV in the model) . Figure 
shows the decay rate for the same mode as a function 
of w, for selected temperatures. The rates become rea- 
sonably insensitive to w above 0.4 meV. The dispersion 
due to the sensitivity on w is a factor contributing to the 
uncertainty of the calculated values. 



V. AMORPHOUS SILICON 

We now present the calculated decay rates for the 1000- 
atom model of a-Si. Due to the computational power lim- 
itations we have sampled the spectrum uniformly with 
about 200 modes for which we computed 2r. The cal- 
culated decay rates are presented as a function of the 
mode frequency for two different temperatures: 10 K in 
Fig. 13 and 300 K in Fig. [S] For comparison the previous 
calculations 6 on a 216-atom model of a-Si are included. 
Overall, the decay rates for the two models agree. The 
rates are on the order of meV (picosecond lifetimes) . Per- 
turbation theory is thus valid for all the sampled modes 
with the exception of few in the lowest part of the spec- 
trum at 300 K (see below). As was shown in Ref. y the 



decay rates as a function of frequency at 10 K (and at low 
temperatures, generally) follow the joint density of states 
Efc/ ^( w 0) — ^(k) — oj(1))] which counts, for a chosen 
mode j, the number of combination decay possibilities 
j — > k + 1 with the constraint of energy conservation. At 
larger temperatures one must add the number of differ- 
ence decay channels j — ► k— I to reproduce, qualitatively, 
the calculated 2T(uj). More detailed physics of the an- 
harmonic decay in glasses and especially the statistics of 
the decay matrix elements can be found in Ref. |6J. 

There are several features which make the calculated 
decay rates for the 216-atom and 1000-atom models 
somewhat different. The first is the overall reduction 
in noise for the 1000-atom model (the data are not 
smoothed as was done in Ref. 0). The reason is both 
the greater model size (spectral averaging) and greater 
w (rectangle averaging). Note that the observed noise 
in the spectrally resolved 2r is about 10% or less, con- 
sistent with a dispersion of the decay rates with w, dis- 
cussed in the previous section. Second, the calculated 
rates for the 1000-atom model are somewhat smaller than 
those of the 216-atom model, that is, the latter model 
appears to be slightly more anharmonic. This is at vari- 
ance with the calculation of thermal expansion^ where 
the 216-atom model seems less anharmonic. The latter 
difference probably can be explained by the anomalously 
large negative mode Gruneisen parameters of the low fre- 
quency resonance modes of the 1000 atom model, as the 
thermodynamic Gruneisen parameter depends on an av- 
erage mode Gruneisen parameter at high temperatures. 
We note that the structural models differ in other ways: 
the smaller model is more topologically constrained^ 
has smaller energy /atom, and has higher density than 
the 1000-atom model. Third, the calculated rates of 
the 1000-atom model extend to a lower frequency region 
as the minimum frequency of the model is smaller than 
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FIG. 9: Calculated decay rates of selected modes in a-Si ver- 
sus temperature. The lines are labeled according to modes' 
frequencies in meV. The lowest frequency mode is a propagon, 
the following two are diffusons (acoustic-like and optic-like) 
and the highest frequency mode is a locon. 

that of the 216-atom model. Finally, some low-frequency 
modes (resonances) at 300 K exhibit giant decay rates, 
comparable to the modes' frequencies. These rates are 
in fact invalid, since they are not consistent with pertur- 
bation theory. However, they indicate what may be ex- 
pected from a full anharmonic calculation (for example 
by molecular dynamics). This important physics issue 
will be discussed elsewhere. 

In Fig. El we plot the temperature dependence of the 
decay rates of selected modes. We show the temperature 
dependence for a propagon, an acoustic-like and an optic- 
like diffuson, and a locon. The low-frequency propagon 
has a divergent lifetime (decay rate vanishes) as temper- 
ature decreases to zero, since there are no two modes 
into which it could decay, due to the energy conservation 
constraint and the existence of the minimum-frequency 
mode. All the other modes have constant decay rates 
at small temperatures. The constant goes smoothly to 
a linear function at large temperatures, which is due to 
the fact that the population density of thermal phonons 
increases linearly with temperature. 



VI. PARACRYSTALLINE SILICON 

The results for the 1000-atom model of p-Si are shown 
in Figs. 1101 and ITT1 which plot 2T as a function of mode 
frequency. For comparison we also present the data for a- 
Si discussed in the previous section. The results are quan- 
titatively similar for both models. There are no anoma- 
lous decay rates appearing in the spectrum of p-Si which 
would be due to the crystalline cluster. In addition to 
the sampling modes, we have computed the decay rates 
specifically for three modes in the mobility edge region 
with the weight at the crystalline cluster greater than 
30%. They are presented in Figs. 1101 and ITT1 by empty 
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FIG. 10: Calculated decay rates of the 1000-atom model of p- 
Si at 10 K. For comparison, the rates for the 1000-atom model 
of a-Si are also shown (shaded area). The empty circles are 
for three modes in the mobility region with more than 30% 
weight on the crystalline cluster, while the inset shows the 
decay rates (filled circles) of modes around 30 meV, which 
have large affinity (weight up to 30 %) for the cluster. 



circles. The decay rates of these modes have the same 
magnitude as those of the other locons. Finally, in the 
insets of the two figures we show the decay rates of modes 
with frequencies around 30 meV, the region of especially 
high affinity for the cluster (see Fig. 0}. Decay rates of 
more than 80 modes in that spectral region are plotted. 
Although many of the modes have large weight (some of 
them up to 30%) on the cluster, most are unbiased. The 
fact that 2T of all of these modes are similar in magni- 
tude at different temperatures implies no special decay 
behavior for the modes of strong affinity for the cluster. 

Figure IT21 shows the temperature dependence of three 
modes with more than 30% weight on the crystalline 
cluster. The modes have frequencies (weight, p) 70,67 
meV (31%,160), 72.68 meV (55%,12), and 73.05 meV 
(59%, 13). In addition, the figure plots the decay rate of 
a "normal" locon with u> = 77.76 meV (0.02%, 8), residing 
outside of the cluster. No unusual features are observed 
in the temperature dependence of the decay rates of these 
modes. Finally, in Fig. 1131 we plot the anharmonic ma- 
trix elements V(J, k, I) of the combination decay j — > k + l 
for the maximally localized mode on the cluster, with fre- 
quency 73.05 meV to visualize the decay channels. The 
figure shows that the dominant channel is a decay into 
two diffusons. Decay into a propagon and a diffuson (the 
points in Fig. 1131 below 15 meV and above 58 meV) is 
somewhat less important. The corresponding matrix el- 
ements are much smaller. This may be related to the 
fact that propagon's weight on the cluster is systemati- 
cally lower than 8.5% (see Fig. [3J- The diffusons' weight 
on the cluster is much more more scattered, with a sig- 
nificant number of diffusons having the weight of 8.5% 
and more. Decay into another locon and a propagon is 




10 20 30 40 50 60 70 80 90 
PHONON ENERGY (meV) 



FIG. 11: Calculated decay rates of the 1000-atom model of 
p-Si at 300 K. For comparison, the rates for the 1000-atom 
model of a-Si are also shown (shaded area). The empty cir- 
cles are for three modes with more than 30% weight on the 
crystalline cluster, and the inset plots the decay rates (filled 
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FIG. 13: Calculated matrix element V(j,k,l) of p-Si for mode 
j maximally localized (59%) on the crystalline cluster [ui(j) = 
73.05 meV] as a function of uj(k). Shown are data for k and 
I obeying energy conservation: \ui(j) — ui{k) — < w/2, 

where w = 0.2 meV (taken to be smaller than w = 1 meV 
used in the calculation, to get a manageable graphics size). 
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FIG. 12: Calculated decay rates for selected locons in p-Si. 
The curves are labeled according to frequency in meV. The 
numbers in the brackets show the modes' weight on the crys- 
talline cluster. 



forbidden by energy conservation. 



VII. CONCLUSIONS 

We have calculated anharmonic decay rates of 1000- 
atom models of a-Si and p-Si using perturbation theory 
with cubic anharmonicity in the interatomic potential. 
The results for a-Si are in agreement with the previous 
perturbative calculations on a smaller model, as well as 
with a molecular dynamics calculation. The results reit- 
erate the previous findings that the vibrational lifetimes 



are on the picosecond time scales, generally increasing 
with increasing frequency. The decay rates of locons are 
idiosyncratic, but are by no means inhibited. Calculated 
decay rates of p-Si are similar to those of a-Si, showing 
little sensitivity to structural properties. These findings 
disagree with the interpretation of recent experiments 
which find decay rates on the order of nanoseconds and 
somewhat greater sensitivity to structural properties. 

The explanation that we offer to account for these dis- 
crepancies is that the calculation and experiment refer 
to two different things. First, as we have pointed out 
earlier, simple (and at such scale usually over-relaxed) 
models like a continuous random network type WWW 
model or a similar model containing a crystalline grain 
used in this study can not faithfully reproduce a broad 
range of various topological features — some or combi- 
nations of which may be responsible for increased decay 
times observed in the experiment — present in a real ma- 
terial. Second, in our calculations only "perturbative" 
decay rates, where a small (infinitesimal) population of a 
single mode goes out of equilibrium are computed. The 
experiments measure the decay of vibrational states ex- 
cited over a large portion of the spectrum. Furthermore, 
the laser excitation produces phonon populations too far 
off the equilibrium to be called small perturbations. In 
Ref. 3, for example, the excited phonon population n lies 
between 0.03 and 0.3, which, for a mode with frequency, 
say, 50 meV corresponds to an effective temperature of 
160 and 400 K, respectively. This is huge compared to 2 
K at which the samples are held. A numerical investiga- 
tion of BickhamA indeed shows that a strong perturbation 
of the vibrational spectrum of a-Si can relax on a 100 ps 
time scale, compared to 10 ps for a weak perturbation. 
In addition to pure vibrational relaxation, it is also likely, 
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as suggested by Bickham and Feldman^ that correspond- 
ingly large local deviations in the atomic displacements 
cause local structural rearrangements which may relax 
to local metastable minima while emitting phonons. It is 
possible that these rearrangements occur on nanosecond 
timescales. 

We conclude that current simple models for a-Si in 
combination with presented above methods of analysis do 
not provide an answer to the question why vibrations in 
real material under the specific experimental conditions 
appear to decay over such long time scales. A better an- 
swer may be obtained by studying larger and more real- 
istic models when computational power required to make 
such calculations tractable becomes available. However, 
considering the consistent results obtained thus far with 



different models (of sizes from 216 to 4096 atoms) and 
different techniques (perturbation theory and molecular 
dynamics) it is rather unlikely that the interpretation of 
the experimental findings in terms of the (pure) vibra- 
tional lifetimes (that is, without considering structural 
relaxation and possibly photoexcited-electron recombina- 
tion processes) will be validated by investigating larger 
models. 
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